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Abstract 

We show that in an SU(2) g) U(l) model with a DSF-like invisible axion it is possible to obtain 
(i) the convergence of the three gauge coupling constants at an energy scale near the Peccei-Quinn 
scale; (ii) the correct value for sin 2 9w(Mz); (hi) the stabilization of the proton by the cyclic Zi^Z^ 
symmetries which also stabilize the axion as a solution to the strong CP problem. Concerning the 
convergence of the three coupling constants and the prediction of the weak mixing angle at the 
Z-peak, this model is as good as the minimal supersymmetric standard model with /isusY = Mz- 
We also consider the standard model with six and seven Higgs doublets. The main calculations 
were done in the 1-loop approximation but we briefly consider the 2-loop contributions. 
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The convergence of the three gauge coupling constants #3, #2, 9i = \/5/3^', and the 
prediction of the electroweak mixing angle are some of the motivations for grand unified 
theories (GUTs) ^, 0, 0, 0. Unfortunately the simplest and more elegant GUTs which 
break into a simple step to the standard model, like SU(5) and some of the 5*0(10) 

and E e GUTs jjJQ, were ruled out by two experimental results. The first one is concerned 
with the fact that using the value of the electroweak mixing angle measured at the Z-peak by 
LEP, i.e., sin 2 8 w{M z ) = 0.23113(15) [a], the running coupling constants extrapolated from 
the values measured at low energies do not meet at a single point or, inversely, assuming 
the convergence of the three coupling constants, the prediction of the value of sin 2 9w(Mz) 
(or, alternatively, we can obtain the convergence point of oti and cti and then the value 
of a s (Mz) can be predicted) does not agree with the experimental value. On the other 
hand, in the minimal supersymmetric extension of the standard model (MSSM for short) 
the coupling constants do intercept at a single point, with a remarkable precision, and the 
predicted value of the weak mixing angle is also in agreement with the observed value f|. 
In fact, this has become one of the most important reason for believing on the existence 
of supersymmetry at the TeV scale and, in particular, on the MSSM. In particular, a new 
scale for physics, /xsusy, is then required. The second trouble with the simplest GUTs is 



'etime. Recent data 
7(. However, this is 



related with the non-observation of the proton decay at the predicted 
by Super-Kamiokande on p — > K + v imply that Tp > 10 33 years P, 
also a trouble for SUSY SU(5) since this theory has d = 5 effective operators induced by 
colored-Higgs triplet that produce a rapid proton decay and it is necessary to appeal to 
fermion mixing in order to keep (tight) agreement with data [9]. Thus, it appears natural 
to ask ourselves if there are options for SUSY yielding the convergence of the couplings, the 
observed value of the weak mixing angle at the Z-pole and an appropriately stable proton. 
The importance of the Higgs boson contributions to the convergence of the couplings has 
been emphasized recently |H| , although exotic fermions can also lead to the gauge coupling 
unification even at the TeV range |Uj. In this cases, unification at lower energy scale is 
possible even without supersymmetry and the proton stability is guar anteed by additional 



assumptions as extra dimensions or dynamical symmetry breaking [ll| , by the conservation 
of the baryon number in the gauge interactions as in the [SU(3)] 3 trinification [10], or m 
[<S77(3)] 4 quartification where the proton decay is mediated only by Higgs scalars |l2j ]. 
Here, we will show that in a multi-Higgs extension of the standard model we have the 



2 



convergence of the three gauge coupling constants at an energy scales of the order of 10 13 
GeV, and the weak mixing angle in agreement with the measured value at the Z-pole. The 
proton decay occurs only throughout dimension 8,9,10 and 11 effective operators because 
discrete Z 13 £*D Z% symmetries forbid all d — 6,7 operators. This model was proposed inde- 
pendently of the issue of the gauge unification, and it has a DFS-like invisible axion 
stabilized against semiclassical gravitational effect by those discrete symmetries [l^j]. The 
use of discrete gauge symmetries in the proton decay problem has been used recently in 
Ref. jl5| in a model with a Z e symmetry, which baryon number violation low dimension 
dangerous effective operators are all forbidden. For other recent uses of symmetries like 
these see ^6|. We compare our results with the usual MSSM showing that the convergence 
of the coupling constants and the prediction of the sin 2 § w (Mz) is, in this model, as good 
as in the MSSM with /^susy = M z . We also consider the SM with seven Higgs doublets, 
and briefly comment the SM with six Higgs doublets and /isusy = 1 TeV. 

In Ref. the representation content of the standard model was augmented, by adding 
scalar fields and right-handed neutrinos, in such a way that the discrete Z13 <S> Z% symmetries 
could be implemented in the model. The particle content of the model is the following: 
Ql = ~ (2,1/3), L L = [yV) T L ~ (2,-1) denote any quark and lepton doublet; 

ur ~ (1,4/3), cIr ~ (1,-2/3), Ir ~ (1,-2), ur ~ (1,0) are the respective right-handed 
components. It was also assumed that each charge sector gain mass from a different scalar 
doublet, hence we have the following Higgs multiplets: four doublets $ u , $d, and $ u [all 
of them of the form (2, +1) = (<p + , <f°) T under SU(2) <S> U(l)} which generate Dirac masses 
for u- and d-like quarks, charged leptons and neutrinos, respectively; a neutral complex 
singlet ~ (1,0), a singly charged singlet h + ~ (l,+2) and, finally, a non-hermitian 
triplet T ~ (3, +2). With the discrete symmetries flavor changing neutral currents are also 
naturally suppressed at the tree level. 

Next, let us consider the evolution equations of the three gauge coupling constants, at 
the 1-loop level, 

a^\M) = a~\M z ) jj cos 2 B W {M Z ) + ^ In (^j , 

a 2 \M) = a-\M z ) sin 2 9 W {M Z ) + ^ In , 

a 3 \M) = ^\M z ) + ^-lnf^-Y (1) 
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where foj are the 1-loop beta-function coefficients 

bi = ~ £ T m {F) + T ^ S ) ~ Y° 2{V) - (2) 

fermions scalars 

For SU(N), we have that Tr = C 2 = iV, with N > 2, for fields in the adjoint representation, 
and Tr(S,F) = 1/2 for fields in the fundamental or antifundamental representation; for 
17(1), C 2 {V) = and T(S a , F a ) = (3/5)Tr(F a 2 /4) for a unification in 577(5) [it is the same for 
the case of SO (10)]. Eq. (J2J) is valid for Weyl spinors and complex scalar fields. Considering 
the extension of the standard model with N g matter generations, Ng Higgs doublets and 
Nt non-hermitian scalar triplets, all of them considered relatively light, 

4 13 

3 9 10 5 

4 1 2 22 

6 2 = - 3 N, + -iV H + 5 iV T - T , 

f>3 = jiV, - 11. (3) 

Only the scalar singlet h + will be considered with mass of the order of the unification scale. 
The evolution equations in Eq. (JTJ implies the unification condition l (Mu) = a\(Mu) = 
«3 1 (Mu) = a^ 1 , which also defines the mass scale M\j\ 



Mu = M z exp 



2tt 



a-\M z ) - faj- 1 (M 2 
|6i + b 2 - §&3 



(4) 



From Eqs. © with = 3, Nh = 1 and iVy = 0, i.e., the standard model, give 
(61,62,63) = (41/10,-19/6,-7). In this case, the evolution of gi is shown in Fig. [TJ In 
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FIG. 1: Running couplings in the standard model. 

this figure (and below) we have used the inputs M z = 91.1876 GeV; a(M z ) = 1/128; and 
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otz{Mz) = 0.1172 [5j]. It is clear that with only the representation content of the SM, there 
is no convergence of the three at a given point Q|. 

On the other hand, with N g = 3,Nh = 4 and = 1 i.e., the model of Ref. jl^ . 
with Eqs. (J3j) we have (&i,&2 ; &3) = (5, —2, —7) and the evolution of the coupling constants 
in this case is shown in Fig. |21 We obtain that the three forces unify at the energy scale 
M v ~ 2.8 x 10 13 GeV and o^ 1 ~ 38. The value predicted for sin 2 9 W (M Z ), using the 




FIG. 2: Running couplings in the present model. 



unification scale as an input, is 



sin 2 # 



w 



-a( M z)(h-b 2 ) In 



and we obtain in this model sin 2 9w(Mz) = 0.2311, which coincides with the measured 
value as it should be since the unification occurs with a very good precision. In order to 
compare this result with others we write at the 1-loop approximation, by eliminating 
the unification scale from Eq. (J2J), 

~ _ ll + lN H + 2N T 



22-\{N H + N T ) 
5 [sin 2 MM z ) 
3 



a(M z ) 



1 



■sin 2 9 W {M Z ) 



(6) 



where b = (63 — 62) /(&2 — b\) llj. For being more general we let Nh and Nt arbitrary to 
see if there are other values for them which could fit with unification. The theoretical ratio 
b defined in the first line of Eq. (jSJ, and which depends mainly on the scalar representation 
content (notice that, at the 1-loop level, b does not depend on N g ), should coincide with 
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the quantity denned in the second line which depends only on the experimental values of 
the coupling constants a and and sin 2 8\y at the Z-peak. The experimental inputs then 
implies b = 0.714 using the second line of Eq. (JBJ). Using the first line of Eq. (JJJ) the minimal 
standard model implies b = 115/218 ~ 0.527 (including the scalar contributions), so that 
it does not match and sin 2 6 'w(Mz) = 0.204 according to Eq. (JHJ) and for this reason the 
unification of this model in non-supersymmetric SU (5) was ruled out by LEP data. In the 
present model we obtain b = 5/7 ~ 0.714 and this value matches in Eq. © and gives, then, 
the observed value for sin 2 Q w {Mz) m as we pointed out above. 

In the MSSM when /zsusy = Mz we have (61,62,^3) = (33/5,1,-3) and the respective 
evolution is shown in Fig. |HJ with Mu ~ 2.1 x 10 16 GeV and the inverse of the coupling 
constant at the unification scale, denoted in this case by a$, is a^ 1 ~ 24. The weak mixing 
angle has also the correct value [l7[ and the same value for b is obtained like in the present 
model. The case when the SUSY scale is of the order of Mz is better than the case when that 
scale is of the order of 1 TeV but we do not show the latter case. As can be seen from Figs. 
2 and 3 and from the value of b, concerning the unification and and the prediction of the 
weak mixing angle, the model of Ref. JjJ is as good as the MSSM and for this reason it was 
not necessary in the present work to take into account of the theoretical and experimental 
uncertainties. 
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FIG. 3: Running couplings in the MSSM with hsvsy = Mz- 



Finally, in the case when N g = 3 but N H = 7 and N T = we have (61,62)^3) = 
(47/10, -13/6, -7), b = 145/206 ~ 0.704 and sm 2 § w {M z ) = 0.230 according to Eq. © and 
the evolution is shown in Fig. HJwith M\j ~ 5.8 x 10 13 GeV j3j. We have also studied the 
case of the SM with six Higgs doublets. In this case the unification of the coupling coincides 
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with that of Ref. [10| and sin 2 §w(Mz) = 0.226. Moreover, the convergence of the three 
coupling constants in the seven Higgs scalar doublets is better than the case of six of such 
doublets Lfl. 
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FIG. 4: Case of the standard model with seven scalar doublets. 

We see that , in order to have the unification of the couplings and the correct value for the 
weak mixing angle at the Z-pole it is not necessary to have low energy supersymmetry. On 
the other hand, four Higgs doublets and a non-hermitian Higgs triplet are more precise, for 
the unification of the coupling constant, than just seven Higgs doublets. As we said before, 
only the singlet h + has been considered having a mass of the order of the unification scale. 
Assuming that this singlet is light we obtain the same value for My but sin 2 Q W {M Z ) = 0.229. 

The contribution of just one Higgs doublet is almost negligible compared to that of 
quarks and leptons when dealing with the renormalization group equations. But in models 
like the multi-Higgs extensions we have discussed here, the Higgs multiplets have a total 
contribution already important to unify the theories. However, we stress again, looking at 
the Eq. (JBJ), that the value of sin 2 9 w(Mz) depends strongly on Nh and Nt selecting in 
this way only few possibilities. Moreover, the prediction of the weak mixing angle in this 
model is not an accident. To see that, with the experimental inputs in the second line of 
Eq. © we have Nh = 7.324 — 3.333iVr and it is clear that the best solution for Nh and 
Nt integers is when Nh = 4 for Nt = 1. The seven doublets model Nh = 7 for Nt = 
is the second best solution. The denominator in the exponent of Eq. (J3J) can be written 
as (l/3)iVfl- + (5/3)A^r + 22, and we see that a larger number of Higgs multiplets imply a 
smaller value for My, however, we obtain the best match condition for b with Nh = 4 and 



Due to the precise measurements in the electroweak sector higher order corrections to 
the 1-loop calculations are important to verify if the solutions above can be stable against 
these corrections. Considering only the gauge coupling constants, their running are 
now solutions of the corrected renormalization group equations 

dai(n) 1 
^ dfj, = 2^ 

The general form of the coefficients bij is given in Refs. For the cases discussed in this 
paper, i.e., the only nontrivial scalar representations under SU(2) are doublets and triplets, 



1 3 

3. Q. 



^f- (7) 



they have the following form 



/ 19 jy I A jV„ _|_ M jy 1 N + 9_ Nm + 22 n 44 jy N 

\N Q + ±N H + 2 4N T fN q + ^N H + ^N T -^ AN q 



(8) 



5^9 1 10^'-" 1 5 1 3 9 1 3 n 1 3 1 3 9 

\ %N g \N g fiV 5 -102y 

The quark top Yukawa coupling being of order of unity, is the only one comparable with 
the OLi. However its contribution to Eqs. (|7|). are unimportant when compared with the 
other contributions in Eq. (jHJ), and does not affect the 2-loops running significantly for the 
values of Nh and considered here. A complete treatment includes, of course, all scalar 
interactions (which in the present model include trilinear interactions) and it is much more 
complicated. However, as an illustration, we will consider only the corrections of the gauge 
coupling constants in Eq. (jHJ). In this case, the numerical solutions to the system of equations 
in Eqs. (J7J) can be found using the Eq. (JTJ) after making the simply substitution bi — » Aj, 
with the Aj extracted numerically. The values of Aj and the respective values for M\j and 
sin 2 §w(Mz), at the 2-loop order, are given in TableE We see from the table that in the cases 
with Nh = 4, N? = 1 and Nh = 7, Nh = 0, the values of the weak mixing angle at the Z 
pole are a little above the experimental value (since this is only an illustration that consider 
just the evolution of the gauge coupling constants it is not necessary to take into account 
the experimental and theoretical uncertainties). This does not rule out these models since 
the correct value of sm 2 §w(Mz) may be obtained at the 2-loop level by constraining the 
self-interactions couplings and the masses of the scalars fields. Hence, this partial analysis 
suggests that solutions in the 1-loop approximation might be stable under higher corrections 
since they are not drastically changed when 2-loop corrections are included. 

Next, we come to the question of the proton stability in the present model. We have 
seen that the energy scale of the unification of the coupling constants is of the order of 10 13 
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Mjj 


sin 2 0w(M z ) 
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1 


5.102 


-1.847 


-7.089 


1.5 


0.235 
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4.659 


-2.213 


-7.089 


5.0 


0.231 
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4.762 


-2.035 


-7.089 


3.4 


0.234 



TABLE I: Values for Aj and sin 2 9\v(Mz) at the 2-loop level for multi-Higgs models. Mjj is in 
units of 10 13 GeV. 



GeV, i.e., smaller than the scale of the non-SUSY SU(5) and near, by a 10-1000 factor, 
to the Peccei-Quinn scale. This is, apparently, a disaster from the point of view of the 
nucleon decay. However, it is not because the model accept discrete symmetries that forbid 
potentially dangerous effective operators of d = 6, 7. With the representation content of the 
model it is possible to impose the following Zis symmetry: 

Q — >■ u 5 Q, u R — > uj 3 ur, d R — > uj^dn, L — > u 6 L, 

V R -> U VR, Ir -> LU4 R , $ m -> UJ^^u, ®d -> V^^d, 

$/ -> cj 2 $«, -> -> v» T ^r 1 ^ 

/i + — >■ uj\h + , 



(9) 

with = e 27rifc / 13 ; = 0,1,. ..,6. Moreover, in order to have an automatic PQ symme- 
try [19], it is also necessary to impose a Z 3 with parameters denoted by cj , £>i, and tuf 1 
with $; transforming with cDf 1 :^, with tDi, while all other fields transform trivially 

under Z 3 . For details see Ref. 

We search for effective operators that are SU(2>)c®SU(2)l®U(1)y and Zi 3 (E)Z 3 invariant. 
No grand unified model is assumed here. Effective operators with d = 6 which induce a rapid 
proton decay for M\j < 10 16 GeV, are given by j^J: 



abcd - e ijk£a(3{dR) C ia {u R ) jb (Q L ) C kac {L L )p d , 



^ abc 
Q (2) 



abed 
(3) 



k£a(3(QL)t aa (QL)jf3b(uR) C kc (^)d, 



^abed - e ijk£a(3£ 1P (QL)taa(QL)j(3b(QL) C k lc (L L )pd, 

°abcd = ei jk (re) a(3 ■ {re) lp {Q L ) c iaa {Q L ) mb {Q L ) c klc (L 1 
(5) 



-bed - e ijk(d R ) c ia {u R ) jh {u R ) c kc (l R ) d , 

(6) 



= ^fcMfaMiftML (Z 



tijk{u R ) c ia (d R ) jh (d R ) c kc {v R ) d , 

(d R )f a (Q L ) jab (v R ) c (Qljkpd, 



(10) 



where i, j, k are SU(3) indices; a, (3, 7 and p are SU (2) indices; and a, b, c and d are generation 



indices. From Eq. 
tensor Dirac matrices 



using Fierz transformations, it is possible to obtain all vector and 
20j . All operators in Eq. (fTU|) are forbidden by the Z 13 symmetry in 



Eq. Q; d = 7 operators formed with those of Eq. ()10|) and the singlet <fi (or 0*) are also 
forbidden. Notice that O^f' 7 ' 8 ^ are also forbidden by Z3. However, there are others B — L 
conserving operators allowed by all the symmetries of the model as 

(1) $ f $ (1) 6 4 (3,4) 6 5 (11) 

of d = 8, 10, 11, respectively, that may induce the proton decay, via four fermion interactions, 
after the spontaneous symmetry breaking. Let us write the proton lifetime as 



151 (12) 

where rp = M^a^mp^ with M 5 is the unification scale in the context of the MSSM, 
M 5 ~ 2.1 x 10 16 GeV; is the respective coupling constant at that unification scale with 
a^ 1 = 24; mp is the proton mass; and ol\j is the coupling constant at Mjj in this model 
with oljj X = 38; finally, £ is a factor depending on the effective operator. Although the d = 8 
operator is suppressed by 1/M^, after the spontaneous symmetry breaking, it induces a four 
fermion interaction proportional to £ = v\v u jMy. Since \v\v u \ % (246 GeV) 2 we have that 
|f I < 7.7 x 10" 23 and since M v /M 5 ~ 1.3 x 10~ 3 , a 5 /a v = 38/24 ~ 1.6, in Eq. JT2J there is 
a factor 10 32 with respect to rp. The d = 10 operator is suppressed by M [/ 2 A 4 and it 

induces four fermion interactions like M^ 2 (v^ / A^O^ where A is a mass scale connecting the 
field <p with the four fermion effective operators A may be My (or the PQ scale) or the 
Planck scale. In this case there is a factor |f |~ 2 = (A/i^) 8 in Eq. (|12|). The enhancement 
on the proton lifetime depends on the scales A and v^. Assuming = 10 12 GeV and 
A = Mpi anc k = 10 19 GeV, we have an enhancement factor of 5 x 10 44 with respect rp. If, 
instead of Mpi an k we use A = Mjj but = 10 9 GeV we still obtain an enhancement factor 
1.7 x 10 24 in the proton lifetime. Finally, if v$ = 10 12 GeV and A = Mjj the proton lifetime 
is raised by a factor two with respect to rp. Similar analysis follows for the d = 11 effective 
operators. Hence, this model survive the proton decay problem since with the natural 
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values of the parameters we have that the proton has a lifetime which is compatible with 
the no observation of its decays at the present experimental level. Moreover, notice that the 
d = 5 effective operator My X LL<& v § v is allowed by the Z 13 symmetry but forbidden by Z 3 . 
However the d = 10 operator M^ 1 A~ 5 LL^ u ^ u (f) 5 gives a Majorana mass to the neutrinos 
with a upper limit of 2 eV, obtained when A = and (<& v ) = 246 GeV. 

Summarizing, we have obtained a multi-Higgs extension of the standard model with 
Zis £g> Z 3 symmetries that imply an automatic PQ, B and L symmetries at the tree level. 
The axion is stabilized against semiclassical gravitational effects by those symmetries and 
they also stabilize the nucleon allowing, at the same time, the unification of the three gauge 
coupling constants at an energy near the PQ scale. Last but not least, the correct value of the 
weak mixing angle at the Z-peak is obtained. Although we can always implement a larger 
Zn by adding more matter multiplets, concerning the unification of the coupling constants, a 
larger number of multiplets or higher dimensional representation of SU (2) affect the running 
of the couplings. Only a limited set of representations is allowed in this respect. We should 
mention that an unification scale near the PQ scale is also obtained in an [SU (3)] 4 model 
but this model has no PQ symmetry in its minimal version. The present model cannot be 
super symmetric at low energy (of the order of TeVs), since the fermion superpartners of 
the Higgs scalars would upset the unification of the gauge couplings, however it is possible 
to have supersymemtry if Hsusy ~ My. It would be interesting to search what sort of 
non-SUSY GUT embed this model. 
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